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ABSTRACT 



Represent a time-shared computer system as a group of N terminals, each 
having submission rate A and exponential (p) task durations, with tasks sub- 
mitted to a central (single) processor. There they are serviced in processor- 
sharing or time-sliced mode. It is shown that the R(t), the response time 
conditional on t , the required processing time, becomes approximately normally 
distributed as t increases. Similar results are derived when N increases. 

Variations of the model consider control: an "inside," processor-shared 
queue services at most c tasks, others queueing first-come first-served 
"outside." Other possibilities are described and analyzed. 

1. Introduction 



The abstraction of computer capacity allocation known as processor sharing 
is an attractive simplification of time slicing, sometimes called round-robin 
scheduling. The idea is well known (Kleinrock, 1976): given that j jobs or 

programs are at the execution stage, each receives service equal to one-j_t± °f a 
time unit per time unit. In other words, if the chance that any single job, 
processed alone, finishes in (t,t+h) is yh-fo(h), (exponential-Markov 
service), then the chance that a particular ("tagged") job in the company of 
(j-1) others finishes in (t,t+h) is p(h/j)+o(h) as h 0. Processor 
sharing of the above type tends to be equitable in that it permits short jobs 
access to processing even if they arrive after, and queue with, longer jobs. 

Apparently the first study of delays to arriving and queueing jobs under 
processor sharing was conducted by Coffman, Muntz, and Trotter (1970) . They 
assumed a steady state M/M/1 system with processor sharing, and were able to 
determine properties of the response time, R , given the processing time 
required by the arriving job. Other papers have also appeared. 
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A recent paper by D. Mitra (1981) analyzes response time, R , under the 
assumption of a closed system. Idealize the behavior of a system of N termi- 
nals and a single computer as a classical machine-repair situation: each think- 

ing terminal (failure-prone machine) applies for computer service at rate A , 
and queued or waiting jobs are served at rate p as long as any jobs are 
present* Markov assumptions are made throughout, so X(t) , the number of jobs 
at the service stage, is a birth and death process with transition rates 



X(t) = j + X(t+h) = j + 1 : Xjh + o(h) (1.1) 

X(t+h) = j - 1 : Pjh + o(h) 

*->- X(t+h) = j : 1 - (Aj+Pj)h + o(h) 

and in particular A. = A (N-j ) , y . = p for j 1 , otherwise being zero. Let 
processor sharing govern service eifort allocation. In Mitra (1981) the dis- 
tribution of response time is characterized, and the moments (e.g. mean and 
variance) are found under interesting conditions, such as that the tagged job 
requiring t time units of processing arrives to find j - 1 accompanying jobs; 
the conditional response time, given only processing requirement t , is given 
particular attention. 

In this paper the previous analysis is generalized and extended. We intro- 
duce the idea of processor sharing in an arbitrary birth and death process en- 
vironment, thus allowing quite general terminal-computer interactions to be 
represented. In the process, the meaning of "system state at the moment of tag- 
ged job arrival** is clarified; see also recent work of Lavenberg and Reiser 
(1981). Response time characteristics are computed under the assumption that 
processor-sharing service rates are processor-state-dependent in a more general 
way than that described earlier; this allows for approximate representation of 
oyerhead penalties and also of job scheduling. Other characteristics of tagged 
job response are also studied, e.g. the accumulated processing work, W(t) , 
actually performed on that job by elapsed time t (t < t = required processing 
time) following job introduction; note that W(R) = t , so the first passage of 
W(t) to t is actually the response time. 

Although differential equations may be obtained for transforms of W(t) 
under various initial conditions, and hence, implicitly, for its distribution, 
the results are far from being explicit and informative. However, central limit 
theorems for additive functionals of Markov processes, or for cumulative pro- 
cesses, allow the conclusion that the accumulated work accomplished by fixed 
time x on a "long** job is approximately normally distributed (Gaussian) . This 
fact in turn allows the conclusion that the response time for a "long** job is 
also approximately normally distributed. Additionally, the normal approximation 
may be shown to be valid for our simple model — and probably for others as well — 
when the number of competing terminals becomes large, i.e. under heavy traffic 
conditions. The quality of the normal approximations for finite job lengths is 
currently being assessed by simulation methods. 

In the latter part of this paper we describe queue control in a processor- 
sharing environment. The expedient is to limit the total number of jobs allowed 
simultaneous processor-shared service at an "inside** queue, with any excess in 
"f irst-come, first-served'* status in an "outside" queue. Long jobs are also 
shifted from inside to outside by a sampling mechanism. It is shown that long 
jobs are favored by a small inside span, c (c being the number simultaneously 
processor-shared), while short jobs are favored by large c . 
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2 . Mean Response Tlmc9 



Begin by describing differential equations for the mean response time to be 
experienced by a tagged, particular, arriving job. Other moments satisfy very 
similar equations. 

(a) Conditioning on Required Time and System State. 

Throughout what follows Markovian assumptions are made: service times 

at the computer are independent and exponential (u) . Generalizations to phase- 
type distributions are apparently possible. 

Let R refer to the response time of a newly arrived job, and 



^(t) = E[R|x(0) = j, W(R) = t] , 



( 2 . 1 ) 



the conditional expectation of response time, given that the job is initially in 
the company of j others (arrives to find j - 1 present) and requires "work" 
or processing time equal to t . Let A . h (resp. y.h) for small h identify 
the infinitesimal generator of the accompanying process, so transition rates are 

as in (1. 1) . 



Consider the possible system changes in (0,h), and subsequently; the 

r ( i ) 

following results occur: let y. = u(j-l) — . • 

l J J ; 



m,(t) = h + m.(t - - 9 - 
j J J 



h)[l - (A +p )h] 



+ A.hm. , , (t - • h) + jT.hm. . (t- 

j J+1 J J 3-1 J 



( 2 . 2 ) 



h) + o(h) . 



The term r(j) represents the fraction of time the processor actually spends 
processing when there are j jobs being processed. 

Allowing h 0 one finds the differential equations 



r(1) 



mj'(t) = 1 - (Aj + Uj)m.(t) + Vj m j +1 (t) + U j m j-1 



(t) 



(2.3) 



This is a standard system of linear differential equations with constant co- 
efficients; initial conditions are nu(0) = ® f° r 3 • 

(b) Conditioning on Required Time. 

If one removes the condition that X(0) = j in accordance with the 
stationary distribution appropriate for an arriving job it follows that the 
expected response time is linear in the required processing time, t . This 
holds for quite general birth-and-death process models, and not just for the 
simple machine-repair setup; see Cohen [1979]. Here is the derivation, in 
outline. 

First, observe that the long-run distribution of X(0) , the number of jobs 
present just after the tagged job enters, is 
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CTT J-l x j-l " c 7 T j pr (J) 



J “ It 2 f • . • , N , 



(2.4) 



where c is selected so that the q f s sum to one. Recall that 

\ \ \ J 

0 1 * * ' 1-1 

ir. » tt « J — is the stationary distribution (assumed to exist) of the 

J 0 u 1 n 2 ...y J 

Markov chain X(t) defined by (1.1) with p = pr(j) . This is intuitively 
apparent, but a formal proof can be based either upon an embedded Markov chain 
formulation, or upon the theory of additive functionals of a Markov process; 
see £inlar ((1975), pp. 269-271). The distribution {q ^ } has also been given 

by Kelly, (1979), p. 12. 

Use (2.4) to remove the condition that X(0) = j; put 



m(t) = E X(0) E[R|X(0), W(R) = t] = . 



(2.5) 



Then in terms of the differential equations (2.3); after multiplying through by 
j/r(j) one obtains (m^(t) = 0) 



m ' (t) ' rh) q j + I_ (yVV 0 + Vj+i (t) + Vj-i (c>1 

N i 

" X r(j) q j ‘ 



( 2 . 6 ) 



Thus it follows that the long-run conditional expected response time is linear 
in the processing time requirement: 



eir|w(r) = t) = t 



N 



l 

j=l 



r(j) 



tE 



X(0) 
r (X(0) ) 



(2.7) 



Apparently no such simple form exists for Var[R|w(R) = t] , although 
Mitra (1981) has given a formula for a particular case. It will be seen, 
however, that the above variance is indeed proportional to t if t is large. 



3. Total Work Completed on a Tagged Job in a Fixed Time 

Turn attention now to W(x) , the total work expended by the computer on 
the tagged job by time x after its arrival, given that the tagged job requires 
exactly t time-units of work for completion. ’ If when the job arrives there 
are X(0) = j customers present, then 
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X(0) = j > 1 . 



(3.1) 



W(t) 



X 




r (X(u) ) 
X(u) 



du, 



(a) The Laplace transform of W(x) . 

Here is the derivation of a differential equation for the Laplace 
transform of W(t) : 

*j(s,x;t) = E{e“ sW(x) |x(0) = j, R = t] , for j > 1 (3.2) 

Now argue in a manner analogous to the discussion prior to (2.2) to write a 
backward equation: for 0 < h < t 

<P .(s,x;t) = (e“ Sr(j)h/j )[l-(A +y )h]<{. (s,x-h;t) 

J J J J 

+ (e -sr(j)h/:i )A j h(}. j+1 (s ,x-h; t) (3.3) 

+ (e sr O) h /jj u ^h(j)j_ i (s,x-h;t) + o(h) , 



leading to 




«. + u. + sr(j)/j)*. + Y j+1 * 



(3.4) 



Initial conditions are 



♦j(s,0;t) = E[e“ sW(0) |X(0) = j, R = t] = 1 , t > 0 , (3.5) 

since initially W(0) = 0 , regardless of the job requirements or the initial 

environment. 

(b) A central limit theorem for W(t). 

Examination of (3.1) shows that W(x) involves sums of contributions 
to work accumulated while the system inhabits various states during the period 
(0 ,t). This suggests that, at least for "long" jobs, i.e. such that required 
processing time t , one can anticipate a nearly-Normal distribution for 

W(t) . An appropriate central limit theorem that establishes this for finite 
birth-and-death models can be found in Keilson ((1979), p. 121); call this 
Theorem K . Alternatively, one can make use of the theory of cumulative pro- 
cesses, see Cox ((1962), pp. 99-101); the latter development is adaptable to 
models more general than the simple bir th-and-death process. 
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In order to apply Theorem K redefine the infinitesimal generator (1.1) to 
describe the behavior of the accompaniment , X' (t) , of the tagged customer; note 
that the relevant generator is now 



X'(t) = j X'(t+h) = j + 1 : Xjh + o(h) 



X'(t+h) = j - 1 



-*■ X'(t+h) = j 

for j = 0,1,2,..., N' = N - 1 . Then 



X j+l h + 



pjh + o(h) 



= P 



j+1 



( . 

1 



U +1 J 



h + o(h) 



1 — (A ! + p ! ) h + o (h) 
J J 



(3.6) 



(3.7) 

(3.8) 



W(t) = / f (X' (u) ) du 
0 



_ r T r(X’ (u)+l)du 



/ 



X'(u) + 1 



and theorem K states that 



W(X) - gT 
o/r 



(D) 



N(0,1) ; 



the constants E, and a 2 are such that 



(3.9) 



(3.10) 



N' 

e = l fa)*: 

j=o J 



(3.11) 



= 2 [ f (0) , f (1) , . . . , f (N ' ) ] 


K 




’f(O) 




0 








"i 


£ 


f (1) 

• 




0 "r . 




• 

• 

_f(N’) _ 



(3.12) 



where in the present case the definition of f(-) is implicit in (3.9), and 
£ is the matrix 



£ = y{[i-A+JU 1 - & 



(3.13) 
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^ being the identity, and 



*0 


*rrl * * * 


*N' 




(3.14) 


*0 


7T^ • • • 


"A* 


9 




"o 


7r| • • • 


"A* 







J^rows are steady-state probabilities for the accompaniment, and fa is defined 
as follows: 



A 0,l = X c/ Y ’ A o,o 1 X c/ Y 



A. ... = X!/y , A. . . = p ! /y > A. . = 1 - v./y 
J.J+1 3 J.J-l J J,J 3 



(3.15) 



v . = A ! + u ! : y = max v . 

J J J 3 ^ 



(V 0 ^0 f V N ! 



(c) A central limit theorem for response time, R(t) . 

A graph of W(x) vs. x starts with W(0) = 0 and increases in 
random straight-line segments until W(x) = t. The value of x at which this 
occurs, x(t) = R(t), is the first-passage time to t of the work process 
{W(x), xj>0} , and is the required response time, so 



P{W(x) <t } = P{R( t) >x } . 



(3.16) 



Now invoke the previous theorem (3,10) concerning asymptotic normality of 
and a standard argument of renewal theo ry, c f . Karlin and Taylor ((1979), 
pp. 208-209) to see that if t = £x + / o 2 tx , then, as t 00 , 




from which it follows that 



W(t) 



R(t) - at 
1 Ft 



_ R(t) - t/£ 



(D) 



A 



N(0,1) 



F 



(3.17) 
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4. Heavy Traffic Analysis of the Response Time of a Processor-Shared Job 



This section investigates the problem of delay of a tagged job requiring 
t units of processing time when it is accompanied by many others, i.e. is in a 
heavily loaded system. Restrict attention to the machine repair model in which 
A . ■ A(N-j) and y = y , and omit the effect of r(j), i.e. r(j) = 1 . Let 

J J -1 

there be N terminals and one processor, with A being the expected terminal 
think time (exponentially distributed), y = Ny 1 being the processing rate of 
arriving jobs; A and y 1 are fixed but N is large and the service rate 
scaling by N is required in order that queue size be of order N . 

Now utilize the fact (Iglehart (1965), and Gaver and Lehoczky (1976)) that 
if X(t) is the number of jobs at the processing stage at t then X(t) can 
be approximated by a diffusion process: 



x(t) = N a(t) + J5 Y(t) ' (4.1) 



where a(t) is a deterministic function of time and (y(t)} is, for the 
present model, a particular Ornstein-Uhlenbeck process. It turns out that when 
N « 



= A(l-a(t)) - v' 



(4.2) 



or 



a(«o) = 1 _ iL 



(4.3) 



which is feasible if A > y' , i.e. under heavy traffic conditions. Furthermore 



d-Y (t) = - AY(t)dt + /A (1-a (t) ) + y' dB(t) , 



(4.4) 



{B(t),t 21 0) being the standard Wiener process. In the long run. 



dY(t) = - AY(t)dt + SW dB(t) . (4.5) 

It is in the environment X(t) described by (4.1) that the tagged job enters. 

It encounters competition for processor-shared service, and so its accumulated 
work completed by fixed time x is essentially 



W(t) . / .aHj. . (4.6, 



Apply the approximation and expand to second order terms in N , the number of 
terminals, to find 









8 



du 



/n Y (u)du 



(A. 7) 



W(t) 



~ / 
0 



Na(u) £l 



Y(u) - | 
l(u) /nJ 



/ 

0 



du 
Na (u) 



- / 



0 (Na(u) ) 



For simplicity, and to enable comparisons with previous results, let 

a(v) ■» a(») , so the tagged job arrives in the steady state. Expression (A. 6) 

then says that for the approximation advanced here, 



W(t) - / 

0 



du 

E[X(-)J 



(E[X (») J ) z 



/ 



0 



Y (u)du , 



0 < t < t 



(A. 8) 



so the expected amount of work done on the tagged job is nearly x/E[X(°°)] , 
and the actual distribution of total work done is approximately Gaussian 
(integral of an Ornstein-Uhlenbeck process), where the Gaussian property results 
from the assumption of many accompanying jobs, and not necessarily because the 
tagged job is long. 

Standard calculations applied to (A. 7) show that, as t -*■ » , 

T x 

*£/ Y(u)du] ^ 0 and Vav[J Y(u)du] ^ (2y’/A 2 )i, so the normal approximation to 
0 0 
accumulated work W(t) has the parameters 



lNa(») ] 



o = 



= 2 m' 

T r 



2u 



lNa(®) J 4 X z [Na («) ] : 



(A. 9) 



from which it follows that the parameters of the normal approximation to R(t) 

are 



a = lNa(”)J cs E[x(~) ] 



e 2 



1 _ 2 U 1 _ 2 f N - E[X(«)] . 

A 7 Na(“) “X 7 E [X(°°) ] X E[X (®) J 1 



(A. 10) 



These formulas state that if think (demand) rate X is very large then, since 
E£ X(°°) ] -*■ N , the variance of response time diminishes, while of course expected 
response time increases like N . This is plausible since in extremely heavy 

t h 

traffic all terminals compete, and the tagged job gets a steady (1/N) — of a 
quantum. For smaller X the expected response time drops with E[x(°°)J , but 
response time variance increases. 

The above derivations are informative but not rigorous. Semigroup methods 
of Berman (1979) can be applied to place the results on a mathematically solid 
basis* Numerical assessment of the results is also of interest. 
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5 * Numerical Comparisons 

I 

In this section a brief investigation is reported of the numerical agree- 
ment between the very simple formulas from heavy traffic theory for the param- 
eters of the accumulated work distribution and those of direct Markov-chain 
cumulative process theory origin. 

Parameters of Total Work 
Examples . Let p - 1, N * 25, 50, with A varying. 



N - 25 



Rates(A) : 


0.01 


0.0222... 


0.05 


0.10 


0.15 


£(M.Ch.) 


0.77 


0.51 


0.16 


0.067 


.055 


(Diffus.) 


— 


— 


0.20 


0.067 


0.055 


o 2 (M.Ch.) 


0. 44 


0.75 


0.22 


.0045 


0.00085 


(Diffus.) 


— 


— 


1.28 


.0040 


0.00079 



N *= 50 


Rates (A) : 


0.01 


0.0222.. . 


0.05 


oao 


0.15 


£(M. Ch.) 


0.53 


0.13 


0.033 • 


0.025 


0.023 


(Diffus.) 


— 


0.20 


0.033 


0.025 


0.023 


cr 2 (M.Ch.) 


0.83 


0.33 


0.0011 


0.000081 


0.000026 


(Diffus.) 


— 


6.48 


0.00099 


0.000078 


0.000025 



The diffusion approximation and Markov chain parameters agree remarkably 
well when traffic is heavy (large A), but, as might be feared, diffusion fails 
miserably for small A . 



6. Queue Control by Service Span and Interruption 

In this section we consider queue control. The central processor now has 
finite service span , c , which may be smaller than the number, N , of ter- 
•minals. This means that if there are i < c jobs in service they are served 
as before "inside" at a rate yr(i) , with processor sharing in effect. How- 
ever, if there are more than c jobs simultaneously requesting service, only 
c of them are served simultaneously, and at rate pr(c) , also with processor 
sharing discipline. The others are queued "outside", with "first-come, first- 
served" service discipline. 
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If there are more than c customers requesting service, the customers that 
are in service "inside" experience independent service interruptions at rate v. 
When service is interrupted, each job in service is equally likely to be moved 
to the end of the queue; thereupon the job at the head of the "outside" queue 
immediately enters service. Both the imposition of the limited processor 
sharing, imposed by c <_ N , and the interruption process are intended to control 
queueing by adjusting the relative attention given to short and long jobs. 

Markovian assumptions are made throughout, so that X(t) , the number of 
jobs requesting service at time t , is a birth and death process with transi- 
tion rates given by (1.1). 

(a) An Auxiliary Process. 

Let R be the response time of a newly arrived tagged job. Since the 
tagged job may not be served until completion when it first enters service, it 
is necessary to introduce an auxiliary process {Y(t); t >_ 0} in order to study 

R . 



In brief summary, if there are i customers (including the tagged job) re- 
questing service at time t and the tagged job is in service, then the state of 
Y(t) is (i,0) . If there are i > c customers requesting service at time t 
and the tagged job is not in service but is in the jth position in queue, the 
state of Y(t) is (i,j) . One can now describe the possible changes in Y(t) 
in a time interval of length h ; details will be presented elsewhere. 

Let 

m (i,j) (t) = EIR l Y(0) = (i,j) ’ W(R) = t] ’ (6.1) 



the conditional expected response time, given that the tagged job is initially 
in the company of (i-1) others and either it is being served inside 
(if j = 0) or it is j_th in the outside queue (if j > 0) . 



ra 



(ij) 



Arguments similar to those of Section 2 yield differential equations for 
(t) which can be solved numerically. A closed-form solution is compli- 



cated and uninformative. It is possible to numerically evaluate the mean 
response time, 



m(t) = E{R|W(R) = t] 



J, <I (i,j)"(i,j) (t) ■ 

1 > J 



whe.e , (lij) 



is the initial distribution encountered by the tagged job. 



The mean response time is not generally linear in t for this model. Note 
that if c - N then this model is equivalent to that considered in Section 2 
and hence as shown is Section 2, m(t) i_s linear in t for that special case. 



(b) An Approximation to Expected Response Time. 

A useful approximation to the expected response time for a job requir- 
ing t units of work is obtained by the following argument. Assume that the 
service rate for the tagged job is the same throughout its processing and is 
equal to the rate that it experiences when it first enters the processor. Thus 
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the tagged job requires t* ■ ~r ( ± ) units of processing time if it enters when 

these are i _< c jobs (including the tagged one) requesting processing. If 

c t 

i > c , then t* “ v . If i > c , then the number of service interruptions 

v 

during t is Poisson with rate — . Each time service is interrupted, the 

C -l 

tagged job spends an expected amount of time (i-c) [v + pr(c)] in queue. 
Thus the expected time spent in queue because of service interuptions is 

(i-c)(— )t*[v + yr(c)] ^ . If i > c , the expected initial wait in queue until 

C -I 

the tagged job starts service is (i-c) [v + yr(c)] . The resulting approxi- 
mation to the expected response time is (see (2. A)) 



c . N 

V it . v ct 

= 1=1 qi r(i) i-c+l Qi r(c) 



( 6 . 2 ) 



N 



£ q^i-c) [v + yr(c) ] 



i=c+l 



vt 
r (c) 



+ 1 



Table 1 gives values for the expected response time, m(t) , and the above 
approximation for various values of A, p, v, c, and t for r(j) - 1 and 
N = 25 terminals. The quality of the approximation (6.2) appears to be 
excellent for all cases considered. 



(c) Numerical Implications. 

Aspects of the behavior of m(t) to be noted from the table are as 
follows. If the amount of processing time required, t , is "small", then 
expected response time is minimized when c is maximized (here c = 25) ; that 
is, when there is maximal processor sharing and no outside queue. If t i- s 
"large", then expected response time is minimized when c = 1; that is, when 
the processor is dedicated solely to the job that is being served, and other 
jobs queue outside in turn. Note that increasing the rate of service inter- 
ruptions by changing v can either increase or decrease the expected response 
time, depending upon job time requirements. 

These behavioral aspects also appear by taking derivatives of the approxi- 
mate average response time A . In particular. 






< 0 
= 0 
> 0 



if 

if 

if 



c " y ’ 

t= 7’ 

t>i. 

v 



If r(j) = 1 
increasing in 



j = !>• 

c if 



,N, then A is decreasing in c if 



t > — , and 
P 



t < — ; A is 
1 p 

A is constant in c it t - — ♦ 
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Finally, arguments similar to those in Section 3 will show that the re- 
sponse time is approximately normally distributed when the required work is 
large. Again, details will be provided in later work. 
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APPENDIX 



To show 
large for the 

First : 

Hence , 

Since, 



Further, 

Hence, 

m = 

where A^xAq = 



[2.7) and (3.10) yield the same value for E([R|w(t) = t] for t 



case Uj = pr(j) . 



q(j) = 



CUjUj 

A q x 



. . x A . , 
T-l 



P 1 X 



j-1 



m 



s l rOT 



= c 



00 . A- X ... X A . , 

V _J 0 Jzl 

r(j) u 1 x ... x y 



Mj = pr(j) 



^ (A 0 x ... x A ) 

m = cy ) j ^ 

j=i p i x ••• x ^ 



CP 



I j(A n x 



n y. j=i 
k=i k 



* vp 



n 

i-j+l 






c = 



n y, 

k=l 1 



j=l x ° * 



* Vi /A 

k=j 



( 

CO 

I 1 (Aq X ... X A ) 

j=l 


00 

n p k 

k=j+l 


00 

1 j 
J . .1-1 


fj-i 1 
n a . 
1=1 x l 


f * 

n p k 

k=i+l 


00 

1 (Aq x ... X A ) 

j=l U J 


00 

lk-j J 
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1 

( 1. 

II r — 3 8 
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j-i 
n a. 
i=i X J 
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00 

" ^ 
[k=j J 





1 by convention. 
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Second : 



X x ... x X 

* + <j) " x — (j+l)7T'(0) j = 0,1,... [X x X = 1] 

1*2 x • . • * ^j+i J. u 



where 



ir*(0) = 



n p, 

k=2 



l j(x x ... x X ) n p 

j=l 3 k=j+l 



«- I Jdgy-.’O) - I ^.'(3-D 

j=0 3 X j-1 J 



« X. X . . . X x . . 

£ r(j) ^ — 

J-1 M 2 X -*- Xy j 



Tr'(0) 



. 00 X, x . . . x X . . 

- 7 — TT ' 

M >1 ^2 X X 



( 0 ) 



since 



Pj +1 = pr ( j+1) . Thus 



c “ 



TT 1 (0) 1 



n M v 

k=2 



" I 

*• j-i 1 



x X 



j-1 



n y, 

k=j 




A comparison of £ with m shows that 
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